Pis'ma V ZhETF 



Evolution of edge states in topological superfluids during the quantum 

phase transition 

M.A. Silaev* and G.E. Volovik +* 

* Institute for Physics of Microstructures RAS, 603950 Nizbny Novgorod, Russia 
+ Low Temperature Laboratory, Aalto University, School of Science and Technology, P.O. Box 15100, FI-00076 AALTO, Finland 
* Landau Institute for Theoretical Physics RAS, Kosygina 2, 119334 Moscow, Russia 



o 

(N 

> 
O 



The quantum phase transition between topological and non-topological insulators or between fully gapped 
superfluids/superconductors can occur without closing the gap. We consider the evolution of the Majorana 
edge states on the surface of topological superconductor during transition to the topologically trivial supercon- 
ductor on example of non-interacting Hamiltonian describing the spin-triplet superfluid "^He-B. In conventional 
situation when the gap is nullified at the transition, the spectrum of Majorana fermions shrinks and vanishes 
after the transition to the trivial state. If the topological transition occurs without the gap closing, the Ma- 
jorana fermion spectrum disappears by escaping to ultraviolet, where Green's function approaches zero. This 
demonstrates the close connection between the topological transition without closing the gap and zeroes in 
the Green's function. Similar connection takes place in interacting systems where zeroes may occur due to 
interaction. 
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1. INTRODUCTION 



General properties of fermionic spectrum in con- 
densed matter and particle physics are determined by 
'topology of the ground state (vacuum). The classifica- 
tion schemes based on topology [Il[2l|3l|4j[5l|6l[7] suggest 
in particular the classes of topological insulators and 
fully gapped topological superfluids/superconductors. 
The main signature of such topologically nontrivial 
'vacua with the energy gap in bulk is the existence of 
.zero-energy edge states on the boundary or at the in- 
terface between topologically distinct domains [SI [H]- 
In Refs. El [7] the classification is based on topo- 
logical properties of matrix Green's function, while the 
'Other schemes explore the properties of single particle 
.Hamiltonian and thus are applicable only to systems of 
'free (non-interacting) fermions. As was found in Rcf. 
, |10j . classifications of interacting and non-interacting 
'fermionic systems do not necessarily coincide. This is 
related to zeroes of the Green's function, which accord- 
mg to Ref. [5] contribute to topology alongside with the 
■poles. Due to zeroes the integer topological charge of 
the interacting system can be changed without closing 
the energy gap, and it is suggested that this may lead to 
the occurrence of topological insulators with no fermion 
zero modes on the interface [TTl [12] . 

In principle, the analogous situation with zero in the 
Green's function and quantum phase transition with- 
out closing the energy gap may also occur in the free 
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Fig.l. Phase diagram of topological states of triplet 
superfluid of "^He-B type in equation in the plane 
(/i, 1/m). States on the line 1/m = correspond to the 
Dirac vacua, which Hamiltonian is non-compact. Topo- 
logical charge of the Dirac fermions is intermediate be- 
tween charges of compact '^He-B states. The line fi — 
marks topological quantum phase transition, which oc- 
curs between the weak coupling '^He-B (with fi > 0, 
m > and topological charge Nk ~ 2) and the strong 
coupling ^He-B (with ^ < 0, m > and Nk = 0). This 
transition is topologically equivalent to quantum phase 
transition between Dirac vacua with opposite mass pa- 
rameter M = ±1^1 . The gap in the spectrum becomes 
zero at this transition. The line 1/m = separates the 
states with different asymptotic behavior of the Hamil- 
tonian at infinity: -ff(p) ±r3p^/2m. The transition 
across this line occurs without closing the gap, 
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fermion case. Example is provided by superfluid "^He- 
B, which belongs to the topologically nontrivial class 
of fully gapped superfluids, which possesses Andreev- 
Majorana fermions on the surface [HI [Til [T51 [TBI [T71 [T51 
[191 [20l [2T] . In the phase diagram of the topological su- 
perfluid/superconductor of the '^He-B type in Fig. [T] the 
topological quantum phase transition (TQPT) between 
two states with different topological charges across the 
line 1/to = occurs without closing the gap [21]. In- 
stead, the asymptotic behavior changes at the momen- 
tum infinity, at p ^ oo, where the Hamiltonian diverges 
and thus the Green's function approaches zero. Such 
scenario is impossible in the models with the bounded 
Hamiltonian jlll 112] , which takes place in approxima- 
tion of finite number of crystal bands. We consider 
the evolution of the spectrum of Majorana fermions on 
the surface of a topological superfluid when the system 
crosses the lines of the TQPT with and without closing 
the gap. 



2. SPECTRUM OF EDGE STATES 

The invariant Nk relevant for '^He-B in Fig. [T]is the 
topological invariant protected by symmetry: 



Nk = 



247r2 



tr 



K I H-^dp^HH 



-'dp,H 



(1) 

where K is matrix which commutes or anti-commutes 
with the Hamiltonian. The proper model Hamiltonian 
which has the same topological properties as superflu- 
ids/superconductors of the ^He-B class is the following 



H = 



P 

2m 



fi \T3- CTl(T ■ p , 



(2) 



where and Ui are Pauli matrices of Bogolyubov- 
Nambu spin and nuclear spin correspondingly; the pa- 
rameter c serves as the speed of light for the Dirac 
Hamiltonian obtained in the limit 1/m — Q and further 
will shall use c = 1. The symmetry K, which enters 
the topological invariant Nk in (H}, is represented by 
the T2 matrix, which anti-commutes with the Hamilto- 
nian; it is combination of time reversal and particle- hole 
symmetries in ^He-B. 

Let us consider the Majorana fermions using the sim- 
plest model of the interface between the superfluid and 
the vacuum, in which the Hamiltonian changes abruptly 
at the boundary |18j . The boundary is at z = and at 
z < we have the equation Hij] = Etp with H = Hq+Hi 
where 



(3) 



Hi = Tl{p.j:(7x + PyCTy) , 

and we use the boundary conditions ip{z = 0) 
Without loss of generality we set py 
Px =P± = ^Jp^ - pi and 

Hx = Tip^a^. 



(4) 

= 0. 
so that 



(5) 



Now it is possible to simplify the equation by choosing 
the wave function transformation U = a^U and V = V, 
where i/j = {U,V)'^. The Hamiltonian then transforms 
as 



Hn = 



P 

2m 



Hi 



-T2P±cry 



(6) 
(7) 



Since ay becomes the good quantum number, this rep- 
resentation allows to reduce the general problem from 
4 (g) 4 to 2 2 matrices. 

Let us consider flrst the solutions corresponding to 
ayip = Tp. Then we get the following equation in the 
form 



2m 



fi] - e Pz + ip± 



which yields the relation between U and V 




= 



2m 



and 



p\ 2 I 

— ^ [I = —TO ± v TO^ + — 2TO/i. 

2to 



The solution is a superposition of two modes decaying 
at z < with Im{pzi.2) > where pli 2 = Pi 2 ^ -PI- 

'4^ = A'iPi+ B^2 ■ (8) 

The boundary condition yields the equation 

Pzi - ip± Pz2 - ip± 



e + ^ — LL e + ^ — 11 

2m ^ 2m ' 



(9) 



For p\_ < 2m/i this equation (O has an exact solu- 
tion 

e=p^. (10) 

The wave functions an '4'2iz) forming the bound 

state have the localization lengths determined by equa- 
tion 

■ ' (11) 



Pzl,2 = Z TO =F Jp'\_ + to2 - 2to^ 



The solutions corresponding to ay^p — —ip yield the 
spectrum e — —p_L, and taking into account the Py de- 
pendence one obtains the helical spectrum of Majorana 

fermions with i/Majorana = c{ayPx - CFxPy) [18_^. 



Evolution of edge states in topological supcrHuids 



3 



(a) |i>m (b) m>|x 




Fig. 2. Spectrum of Andreev-Majorana fermions, 
localized states on the surface of topological super- 
fluid/superconductor of the ^He-B class (red solid lines) 
for (a) > m > and (b) m > jj,. The spectrum of 
bound states terminates when it merges with contin- 
uous spectrum in bulk (yellow color), whose border is 
shown by blue dashed line. 

For m > the spectrum of Andreev-Majorana 
fermions e = ±p± is shown by the red solid line in Fig. 
Ofor e > 0. The bound states are confined to the region 
< ■y/2m/i. They disappear when their spectrum 
merges with the continuous spectrum in bulk. The edge 
of continuous spectrum is shown by blue dashed line in 
Fig. [21 li m > p the minimum of the bulk energy 
spectrum increases monotonically with momentum p± , 
therefore the bulk gap is 

A = ^ , m> fi. (12) 

If /i > TO the minimum of the bulk energy is non- 
monotonic function of p± having the minimum at 
pTnm _ ^2m{ijL — to) where the bulk gap is 

A = ^J2^m - to2 , < TO < /i . (13) 

The line m = fj, marks the non-topological quantum 
phase transition - the momentum space analog of the 
Higgs transition, when the Mexican hat potential as 
function of p± emerges for /i > to [6] . 

3. EVOLUTION OF EDGE STATE AT 
TOPOLOGICAL QUANTUM PHASE 
TRANSITION 

Let us first consider the behavior of the spectrum 
of Majorana fermions at the topological transition at 
which TO crosses zero. When to approaches zero, to — > 0, 
the region of momenta where bound states exist shrinks 
and finally for to < 0, i.e. in the topologically trivial 
superfluid, no bound states exist any more. Simultane- 
ously the gap in bulk, which at small m is A w -v/Sro/I 
according to Eq. (|13p . decreases with decreasing m and 



nullifies at to = 0. This corresponds to the conven- 
tional scenario of the topological quantum phase tran- 
sition, when at the phase boundary between the two 
gapped states with different topological numbers the 
gap is closed. The same happens at the TQPT occur- 
ring when fj, crosses zero (see phase diagram in Fig. [T|) . 

Now let us consider what happens with bound states 
in the case if the TQPT occurs in the opposite limit, 
when TO changes sign via infinity, i.e. when 1 /to crosses 
zero. This topological transition occurs without closing 
of the gap. In this case the bound states formally ex- 
ist for all Px even in the limit 1/to — > 0. However, in 
this limit the ultraviolet divergence takes place. Two 
components in the superposition (jS]) for the wave func- 
tion of Majorana fermion have characteristic lengths 
determined by imaginary momentum in Eq. (|lip . At 
1 /to — >■ these lengths become 

L^l = Im p^i w ^ - P^/(2to) /i , 

L^2 ~ I™ Pz2 ~ 2to — > oo . (14) 

The length of the first component remains finite in this 
limit but the dimension of the second component shrinks 
to zero and thus leaves the region of applicability of the 
model Hamiltonian. As a result the wave function of the 
bound state cannot be constructed any more. In other 
words, if the TQPT from topologically non-trivial to 
the trivial insulator (or superconductor) occurs without 
closing the gap, the gapless spectrum of surface states 
disappears by escaping via ultraviolet. 

4. CONCLUSION 

We considered two scenarios of the evolution of the 
spectrum of the gapless edge states at TQPT. One sce- 
nario refers to the traditional case, when the gap is nul- 
lified at the transition. In this case, when the TQPT 
from the topological state to the topologically trivial 
one is approached the spectrum of Majorana fermions 
shrinks to zero and vanishes after the transition. 

The other scenario takes place, when the TQPT oc- 
curs without the gap closing. In this case, the spectrum 
of Majorana fermions vanishes by escaping to ultravio- 
let. The characteristic momentum of one of the wave 
functions relevant for the forming of the bound state di- 
verges at the TQPT as Im pz2 = 2to oo. This limit 
corresponds to formation of zero of the Green's function, 
G — — H) — > 0. Thus, similar to the interacting 

systems, the two phenomena - TQPT without the gap 
closing and zeroes in the Green's function - are closely 
related. That is why we expect that the same scenario 
with escape to the ultraviolet takes place for the inter- 
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acting systems; if due to zeroes in the Green's function 
the TQPT in bulk occurs without closing the gap, the 
spectrum of edge states will nevertheless change at the 
TQPT, and this change occurs via the ultraviolet. 

In future it will be interesting to extend the con- 
sideration to the bulk-vortex correspondence ^31 Ull 
[25l [26l [27] . i.e. to study the evolution of the spec- 
trum of fermion zero modes localized inside the core 
of vortices (or other topological defects) in the process 
of topological quantum phase transition in bulk with 
and without the gap closing. The escape via the ultra- 
violet can be also considered for the gapless topological 
matter with the fermion zero modes forming flat bands 
[281 [291 [3Ql [31] and Fermi arcs [211 [3l [33]. 
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